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We present new magnetic field measurements of the K2 main-sequence star e Eri based on principal components analysis 
(PCA) line-profile reconstructions. The aim of this paper is to quantify the surface-averaged magnetic field and search 
for possible variations. A total of 338 optical echelle spectra from our robotic telescope facility STELLA with a spectral 
resolution of 55,000 were available for analysis. This time-series was used to search for the small line-profile variations 
due to a surface magnetic field with the help of a PCA. Evidence for a spatial and temporal inhomogeneous magnetic field 
distribution is presented. The mean, surface averaged, magnetic field strength was found to be (B) = 186 ± 47 G in good 
agreement with previous Zeeman-broadening measurements. Clear short-term variations of the surface averaged magnetic 
field of up to few tens Gauss were detected together with evidence for a three-year cycle in the surface-averaged magnetic 
field of € Eri. 
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1 Introduction 

eEri is among the bright stars that are mildly magneti¬ 
cally active and were a planet detection had been claimed. 
Low-amplitude radial velocities variations of eEri were 
first reported by Campbell et al. (ll988Lpatzes et al. (2000)| 
had announced a planet around eEri with an orbital pe¬ 
riod of 2500 d and an eccentricity of 0.6. Benedict et al. 
(12006b and Reffert & Quirrenbach (1201 lb refined the or¬ 
bital solution and combined the radial velocities with as¬ 
trometric observations and derived a likely mass for the 
planet of 1.55 Mjup. However, Zechmeister et al. (12013b 
could not confirm the existence of this planet from higher- 
precision radial-velocity data nor did Anglada-Escude & 
Butler (12012b succeed to coherently combine the different 
radial velocity data sets with the original 2500-d period. 

Detailed studies of the basic parameters of eEri were 
derived from high-resolution spectroscopy by Steenbock 
& Holweger (11981b . Drake & Smith (11993b . Valenti et 
al. (11995b . Riiedi et al. (11997b . among others, and from 
optical interferometry by Baines & Armstrong (12012b . 
These studies converged on the picture that eEri is a 
single, moderately-young (ft!800Myr), K2 main-sequence 
star with an effective temperature of 5100K, a mass 
of 0.82 Mq, a radius of 0.74 Rq, subsolar metallicity of 
—0.17, and barely resolvable rotational line broadening 
(usini = 1.6kms“^). The stellar rotation period was first 
determined from the Mt. Wilson H&K data to ll.lOd by 
Gray & Baliunas ( 11995b . its range further discussed in Don¬ 
ahue et al. (11996b . Time-series photometry was obtained by 

* Based on data obtained with the STELLA robotic telescopes in Tener¬ 
ife, an AIP facility jointly operated with lAC. 
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Erey et al. (11991b and from the 12 spots detected with a 
rotation period of between Ril0-12d a differential rota¬ 
tion very much equal to the Sun’s was derived. Croll et al. 
(12006b . Erohlich (12007b . and Savanov ( 12009b analyzed pho¬ 
tometric data from the MOST satellite and found two spots 
or spot groups separated by 180° in longitude and individ¬ 
ual spot rotation periods of 11.35 d and 11.55 d from which 
a differential surface rotation of roughly half the solar value 
was derived. 

Gray & Baliunas (11995b found evidence of a 5-year ac¬ 
tivity cycle from Mt. Wilson H&K observations. The cycle 
period was revised just recently by [Metcalfe et al. (201 3)| 
from a 45-yr long Ca II H&K data set based on Mt. Wilson, 
CTIO, Lowell, CASLEO, CPS, and HARPS observations. 
They found two coexisting 2.95±0.03yr and 12.7±0.3yr 
periods with evidence for a Maunder minimum-like state in 
the late 1980’s followed by the resurgence of a coherent 3- 
year cycle. 

The line broadening of e Eri due to a magnetic field was 
first announced by Marcy (11984b and Gray (11995b . Marcy 
(11984b inferred a surface magnetic field between 620 G and 
2850 G and a filling factor, /, between 0.89 and 0.2, respec¬ 
tively. These values were confined to 1 kG and / = 0.35 
in subsequent work by Basri & Marcy (11988b and Mary 
& Basri (11989b . Gray (11995b obtained a magnetic field 
strength of 1.9kG and / = 0.3. Saar (11988b improved the 
detections by consequently applying radiative-transfer tech¬ 
niques in magnetized atmospheres but still obtained rather 
different magnetic-field values between 0 and 3.5 kG from 
different spectral lines. He, among others, pointed out the 
limitation due to substantial systematic effects in the data, in 
particular due to line blending. Eirst infrared magnetic-field 
measurements from Zeeman broadening were presented by 
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L. T. Lehmann et al.: Magnetic field measurements of e Eridani from Zeeman broadening 



Fig. 1 STELLA example spectrum of eEri. Shown is the full echelle order #92 from 6140-6230 A containing the 
magnetic-sensitive Le I line at 6173 A. The spectral resolution is 55,000. 


Valenti et al. (119951 1 and yielded a photospheric coverage 
of just 8.8% with a 1.44kG field. Riiedi et al. d 19971 1 re¬ 
discussed the impact of blending and confined the spa¬ 
tially averaged surface field strength on eEri to 165±30G 
and concluded that it is difficult (at best) to separate / 
and B. Finally, Kochukhov et al. (1201 11 1 presented mod¬ 
ern linear and circular polarization observations with the 
new HARPS spectropolarimeter. The Stokes-V profiles of 
e Eri changed systematically from night to night throughout 
the rotation period of the star. The mean longitudinal mag¬ 
netic field varied between —5.8 ± 0.1 G and 4-4.7 ± 0.1 G. 
Recently, [Jeffers et al. (2014j| could also detect variations 
in the large-scale magnetic field geometry. They examined 
the evolution of the magnetic field topology of eEri us¬ 
ing spectropolarimetric observations and Zeeman-Doppler 
imaging for six observational epochs spread over a period of 
nearly seven years. They found similar S-index variation for 
the last seven years as [Metcalfe et al. (2013)| including the 
hint of a three year cycle. Using Zeeman-Doppler imaging 
[Jeffers et al. (2014j| showed a highly variable magnetic field 
geometry. The poloidal component varies from strong dipo¬ 
lar to monopolar and the toroidal component varies from 
being non-existent to being the dominant component of the 
magnetic field energy. 

In this paper, we present and analyze time-series opti¬ 
cal spectroscopy of eEri. We employed our robotic tele¬ 
scope STELLA-II and its echelle spectrograph SES in 
Tenerife. Section [2] presents and describes the new data. 
Sect. [3 describes the Zeeman-broadening PCA technique 
and Sect. [4] presents the results. Our conclusions are pre¬ 
sented in Sect. [5] 

2 Observations and data reductions 

Echelle spectroscopy was taken with the 1.2 m STELLA 
telescopes between July 2008 and October 2013. Spectra 
were taken with an exposure time of 300 s and achieved a 
signal-to-noise (S/N) ratio between 400:1 to 100:1 per res¬ 
olution element depending on weather conditions. A total 
of 338 spectra were used in this analysis. STELLA-I and 
STELLA-II are fully robotic telescope that make up the 
STELLA observatory at the Izana ridge on Tenerife in the 
Canary islands (Strassmeier et al. l2004ir2010l) . The fiber-fed 
STELLA Echelle Spectrograph (SES) is the one telescope’s 


only instrument. It is a white-pupil spectrograph with an 
R2 grating with two off-axis collimators, a prism cross dis¬ 
perser and a folded Schmidt camera with an E2V 2kx2k 
CCD as the detector (until mid 2012). All spectra have a 
fixed format on the CCD and cover the wavelength range 
from 388-882 nm with increasing inter-order gaps near the 
red end starting at 734 nm towards 882 nm. The resolving 
power is i?=55,000 corresponding to a spectral resolution 
of 0.12 A at 650 nm (3-pixel sampling). An example spec¬ 
trum is shown in Fig. [T] We note that the SES received a 
major upgrade in summer 2012 with a new cross disperser, 
a new optical refractive camera, and a new 4kx4k CCD. A 
bit earlier, the SES fiber was moved to the prime focus of 
the second STELLA telescope in 2011. Further details of 
the performance of the system were reported by Weber et 
al. (120121) and Granzer et al. (120101 ). 

Data reduction is performed automatically using the 
IRAF0-based STELLA data-reduction pipeline (Weber et 
al. l201 lb . Images were corrected for bad pixels and cosmic- 
ray impacts. Bias levels were removed by subtracting the 
average overscan from each image followed by the subtrac¬ 
tion of the mean of the (already overscan subtracted) master 
bias frame. The target spectra are flat fielded with a nightly 
master flat which has been normalized to unity. The nightly 
master flat itself is constructed from around 50 individual 
flats observed during dusk, dawn, and around midnight. Af¬ 
ter removal of the scattered light, the one-dimensional spec¬ 
tra were extracted with the standard IRAF optimal extrac¬ 
tion routine. The blaze function was then removed from 
the target spectra, followed by a wavelength calibration us¬ 
ing consecutively recorded Th-Ar spectra. Finally, the ex¬ 
tracted spectral orders were continuum normalized by di¬ 
viding them with a flux-normalized synthetic spectrum of 
the same spectral classification as the target in question. 

3 Magnetic field measurement 

3.1 PCA reconstruction of the Zeeman pattern 

The principal component analysis (PCA) or Karhunen- 
Loeve transformation ( [Bishop 1995| ) is a multivariate algo- 

* The Image Reduction and Analysis Facility is hosted by the Na- 
tional Optical Astronomy Observatories in Tucson, Arizona at URL 
iraf.noao.edu. 
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rithm, that is able to analyze tiny variations in time series. 
The PCA algorithm decomposes the entire data set of ob¬ 
served Stokes I profiles into a new coordinate system, which 
is of advantage to detect the directions of high variances (see 
[Martinez Gonzalez et al. (2008)| and [Carroll et al. (2009)l l. 
The new basis has to explain the largest amount of vari¬ 
ances of the data set by using as few as possible orthogonal 
basis vectors. The PCA computes this set of new basis vec¬ 
tors by calculating the eigenvectors of the covariance matrix 
of the data. The eigenvectors, also called eigenprofiles, are 
ordered by their associated eigenvalues. Therefore, the PCA 
will project the systematic and most coherent features to the 
first few eigenprofiles and will separate them from noise and 
from each other. 

With a set of computed eigenprofiles fi(A), one can re¬ 
construct the Stokes I profiles I(f, A) by 


I(f,A) = ^c,(f)f,(A) 


( 1 ) 


where Ci{t) = I{t, A)fi(A) is the scalar product between the 
Stokes I profile I(t, A) and the eigenprofile fj (A). The coef¬ 
ficient Ci{t) is therefore a measure for how strong a certain 
feature, characterized by an eigenprofile fi(A), is present in 
the Stokes I profile given at time t. 

It is well known that the magnetic field affects spectral 
lines in a certain way. An example is shown in Fig. |2] We 
simulate the Fel 6173 A line without magnetic field and 
with a magnetic field of 1300G on conditions similar to 
eEri. It is clearly visible that the magnetic field leads to 
broadened line wings and thus decreased line depth. The 
variations in the activity of e Eri will produce tiny variations 
with similar shape in the spectral lines. To detect such little 
specific variations in the time series we used the PCA. 


Skumanich & Lopez Ariste (2002) demonstrated the 


physical content of PCA applied to Stokes profiles. They 
used the weak field approximation, which assumes that the 
magnetic field is weak enough to represent the Stokes I pro¬ 
files by the low order terms of its Taylor series. 


where 


I{t, A) = Io(f, A) + a{t)dIoit, X)/dX 
P{t)dll{t,X)/dX^ + ..., 

Io(f,A) 


( 2 ) 

co(f)fo(A). 

Skumanich & Lopez Ariste (2002) show that the coef¬ 


ficient a{t) could be associated with the line-of-sight 
velocity uios and /3(f) with the magnetic Zeeman splitting 
AAz. With the help of [Jefferies et al. (1989)[ they showed 
that 


/3(f) = (1 -I- cos^ (/))AA|/2, 


(3) 


where AAz = 




rB{t), SO that /3(f) 


_ Skumanich & Lopez Ariste (2002)[ and 

[Jefferies et al. (198^ emphasized that the variations 
caused by magnetic fields will have the shape of the second 
derivative of the Stokes I profiles and could be analyzed by 
its coefficients. 



Fig. 2 Synthetic line profiles of the magnetic-sensitive 
Lei line at 6173A without a magnetic field (dashed line) 
and with a surface averaged magnetic field of 1.3 kG (full 
line). 


To prove this prediction, we create a set of 101 synthetic 
Le I 6173 A line profiles which vary in line-of-sight velocity 
from0-500ms“^ and in the magnetic field from0-300G. 
In Lig. the average mean profile and its first and second 
derivative is shown, Lig. (the row below) shows the first 
three eigenprofiles computed by the PCA. The first and sec¬ 
ond eigenprofile matches very well with the first and second 
derivative, as the theory predicted. The variations in veloc¬ 
ity and magnetic field are orthogonal to each other and can 
be separated by PCA. Additionally, the zeroth eigenprofile 
shows the shape of the average synthetic line, as expected, 
because the greatest variation of the Stokes I profiles to a 
zero-line profile is the average line profile. Therefore, the 
PCA algorithm enables us to detect tiny variations caused 
by the magnetic field in our time series of 338 Stokes I spec¬ 
tra. 

3.2 Selection of spectral lines and denoising 

Spectral lines are selected from the VALD-II database 
(Kupka et al. 119991 1. The selection focused on lines of ap¬ 
proximately equal strength, as unblended as possible, and 
with a high Lande factor. As threshold we set 0.8 in the line 
depth and a minimum Lande factor of 1.59. A total of 30 
such lines were selected and are listed in Table [T] Lor these 
lines we expect a comparable dependency on other line¬ 
broadening mechanisms. All 30 individual spectral lines 
were transformed into velocity space and averaged to ob¬ 
tain mean spectral lines profiles. Therefore, the S/N-ratio 
increased by approximately a factor of five. In Lig. 01 we 
show the 338 normalized spectral lines profiles we used as 
our Stokes I profiles for analysis. 
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Fig. 3 Synthetic line profiles analysis for Fel 6173 by 
varying the line-of-sight velocity from 0-500 m s“^ and the 
magnetic field from 0 - 300 G. a From left to right the av¬ 
erage synthetic profile and its first and second derivative, b 
The first three eigenprofiles computed with PCA. 


Table 1 Spectral lines used in this analysis. 


El. 

A (A) 

Lande 

factor 

El. 

A (A) 

Lande 

factor 

Cri 

4545.953 

1.92 

Eel 

4903.310 

2.24 

Cri 

4591.391 

2.00 

Eell 

4923.927 

1.69 

Eel 

4602.941 

1.74 

Eel 

4938.814 

2.01 

Eel 

4611.279 

1.89 

Cri 

4942.496 

1.85 

Cri 

4616.124 

1.67 

Cri 

4964.927 

2.17 

Cri 

4626.173 

2.01 

Eel 

4967.897 

1.88 

Eel 

4635.846 

2.09 

Eel 

5044.211 

1.77 

Eel 

4673.270 

1.84 

Cri 

5348.315 

1.59 

Eel 

4690.138 

2.16 

Mnl 

5394.677 

1.86 

Eel 

4704.948 

2.49 

Mnl 

5432.546 

2.14 

Mnl 

4754.042 

1.65 

Eel 

5497.516 

2.26 

Eel 

4768.396 

1.62 

Eel 

5501.465 

1.88 

Mnl 

4783.427 

1.97 

Eel 

5506.779 

2.00 

Eel 

4848.883 

2.01 

Eel 

6213.430 

2.00 

Eel 

4878.211 

3.00 

Eel 

6219.281 

1.66 


3.3 Eigenprofiles analysis 


To make the PCA more sensitive to small variations in the 
wings of the spectral lines profiles, we analyzed the inverted 
(1— Stokes I) profiles. 


Fig. |5] shows the comparison of the theoretical expec¬ 
tations (see Sec. 13.Il l with our analysis results. The first 
row. Fig. EJi, shows the average profile over all spectra 
in the left panel, and its first and second derivative in 
the other panels. It illustrates the pattern we are look¬ 
ing for. The second row. Fig. |3), shows the zeroth, first, 
and second eigenprofile of the spectral lines profiles de¬ 
rived with PCA. The zeroth eigenprofile in the left col¬ 
umn represents the mean profile of the 338 spectral lines 
profiles. As expected, the variations of the second eigen¬ 
profile match the second derivative of the observed aver¬ 
aged profile. The variations are in good agreement with a 
correlation coefficient of 82 %. This correlation indicates, 
following Skumanich & Lopez Ariste (2002) variations in 



Velocity [km/s] 


Fig. 4 The 338 normalized spectral lines profiles used for 
magnetic-field analysis with PCA. 


the magnetic field affecting the spectral lines on eEri. By 
comparison, the first eigenprofile bears resemblance with 
the first derivative. According to Sect. 13. II this eigenprofile 
indicates a velocity drift of the 338 spectral lines profiles 
against each other and were not further examined. The first 
and second eigenprofiles are orthogonal to each other, so 
that the effects from the velocity drift are completely sepa¬ 
rated from signals resulting from magnetic field variations. 
This proves that the variations detected reside mostly in the 
second eigenprofile due to the joint Zeeman pattern. The 
eigenprofiles in higher orders are containing noise. 

We note that in our PCA all available spectral line pro¬ 
files of different wavelengths and times have to enter into 
the analysis simultaneously in order to achieve the disen¬ 
tangling of the individual eigenprofiles into patterns which 
represent the first and the second spectral derivative (and 
thus a velocity and magnetic field variation). We experi¬ 
mented with two subsets of the given set of spectral lines, 
i.e., we separated the very high Lande factor lines (>2.0) 
from the low and moderate Lande factors (<2.0) and per¬ 
formed the same analysis on the two data sets as before. 
Neither of the two data sets could reconstruct an eigenpro¬ 
file which matches the second derivative with an acceptable 
correlation, nor could such an eigenprofile be found within 
the most significant ones. This is not a problem of the new 
technique but of the data quality, the split data simply do not 
allow such a low detection threshold anymore. 

The coefficients of the second eigenprofile C 2 (f), shown 
in Fig. |6| can now be used to quantify the magnetic field 
and its variations. A coefficient of C 2 (f) ~ 0 represents the 
average magnetic field at time t, while positive and negative 
coefficients represent accordingly weaker or stronger fields. 
The coefficients in Fig.|6]show roughly a long-term oscilla¬ 
tion around zero. We also see scattering of the coefficients 
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b. 

0th eigenprofile 



1st eigenprofile 2nd eigenprofile 




Fig. 5 Comparison of theoretically expected profile patterns with results from PCA. a The observed averaged spectral 
lines profile and its first and second derivative, b The first three eigenprofiles of the spectral lines profiles. 


Year 

2009 2010 2011 2012 2013 2014 



Fig. 6 The coefficients of the second eigenprofile of the 
338 spectral lines profiles derived with PCA. 

on small time scales. The single seasons are clearly sep¬ 
arated by different seasonal mean coefficients. In the next 
step, we will convert the coefficients to values of magnetic 
field. 

3.4 Magnetic field determination 

The actual quantification of the magnetic field, is done 
through a least-squares comparison of the coefficients from 
the observations with synthetic profiles of known B • Vi. 
We assumed / = 1 to obtain a surface average magnetic 
field and generated 1,000 synthetic profiles with B • Vi be¬ 
tween OG and 333 G to analyze this series with PCA. The 



Fig. 7 The calibration curve used to convert coefficients 
of the second eigenprofile C 2 [t) to magnetic field values ■ 
f. Every tenth coefficient of the 1,000 synthetic profiles is 
plotted and the corresponding linear fit function. 


resulting C 2 {t) coefficients relate linearly with {B ■ V7)^^ 
see Eq. and a simple regression then defines the trans¬ 
formation to surface-averaged magnetic field strength. With 
the help of the calibration curve, Eig.|7] we convert the co¬ 
efficients of the second eigenprofile C 2 (t) to magnetic field 
values B • Vi- Hereinafter we defined for the surface aver¬ 
age magnetic field 

= (4) 
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Year 

2009 2010 2011 2012 2013 2014 



Fig. 8 Surface averaged magnetic field measures for the observing seasons 2008/09 to 2013/14. The solid line indicates 
a sine-fit using a period of 1096 d (RiS.O yr), derived by a Lomb-Scargle periodogram. 


4 Results 

4.1 The magnetic field of eEri 

Fig. [8] shows all 338 individual {B) measures with its un¬ 
certainties. We measure a mean magnetic field over all sea¬ 
sons of 186±47G, where the error represents the standard 
deviation. The mean magnetic fields of single seasons are 
listed in Table |2] We see that the average magnetic field of 
the single seasons is oscillating with maxima in the seasons 
2009/2010 and 2012/2013 and minima in 2011/2012 and 
2013/2014. The variations in single seasons are largest in 
the seasons 2011/2012 and 2012/2013 and smallest in the 
season 2013/2014. In the very end of season 2009/2010, 
we clearly see a separated group of few observations where 
magnetic field values are higher as in the rest of the season. 
Season 2010/2011 consists of only 10 observations, which 
makes the variation analysis not reliable, but the data points 
fit in the overall trend. 

Another issue from Fig.|8]is that the magnetic field val¬ 
ues vary over few tens of Gauss at significantly shorter 
time scales then the rotation period of eEri. Their aver¬ 
age internal errors are ±13 G. Short-term variations on a 
1 - 5 cr level are evident within an observing season. This 
could be an indication that the topology of eEri’s mag¬ 
netic field is spatially inhomogeneous. The strong variations 
could appear from local active regions which rotate in and 
out of the line of sight. This would be in a good agreement 
to the large-scale magnetic field geometry determined by 
[Jeffers et al. (2014)| A periodogram analysis does not re¬ 
veal a significant short period, which could be interpreted 


Table 2 Seasonal values of magnetic field measurements. 


Season 

Number of 

observations 

Mean magnetic field 
{B) (G) 

2008/2009 

45 

206 ± 15 

2009/2010 

55 

229 ± 23 

2010/2011 

10 

200 ± 13 

2011/2012 

98 

165 ± 32 

2012/2013 

65 

230 ± 21 

2013/2014 

65 

124 ± 25 

all 

338 

186 ± 47 


as an indicator for temporal inhomogeneities of the active 
regions. 


Temperature surface inhomogeneities like starspots may 
contribute to the disk-integrated line strength, and thus its 
broadening, even if the spots remain unresolved. Because 
in the present paper it is our intention (like in most other 
Zeeman-broadening works) to determine the surface aver¬ 
aged magnetic field, one of the built-in limitations is the 
still common (but surely simplified) homogeneous two- 
component photospheric model, i.e. two magnetic compo¬ 
nents described by a single temperature model (see also 
Sect. 6.4 in Valenti et al. (1 19951 1 for a detailed discussion 
of the possible error budget of this assumption). 
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Periodogram 



0.002 0.006 0.010 0.014 

Frequency (cycles/day) 


Fig. 9 Lomb-Scargle periodogram for Fig. |8] The dashed 
line represents a FAP of 10“®. The vertical dashed-dotted 
line indicates the maximum peak. 


4.2 Time variability 


Fig. [U shows a sine-curve fit to the magnetic-field data that 
indicate a period of 3.0yr (solid line). In Fig. |9] the corre¬ 
sponding periodogram is given. 


[Hatzes et al. (2000)| were the first to report a 
3-years magnetic activity cycle in eEri. Recently 
[Metcalfe et al. (2013)| suggested the simultaneous ex¬ 
istence of two cycle periods of 2.95±0.03years and 
12.7±0.3years. They used the Call H&K S-index to 
estimate the magnetic activity cycles. Our time series of 
{B) measurements of eEri show a similar behavior as the 
published data of the Call H&K S-index in Eigure 1 by 
[Metcalfe et al. (2013) and the published data of the S-index 
in Eigure 2 by Jeffers et al. (2014)] 


We analyzed our time series for long-term periods with 
the help of a Lomb-Scargle (Scargle 1982) analysis and 
determine a period of 1096±471d (3.0±1.3yr). To esti¬ 
mate its false alarm probability (EAP) we used two dif¬ 
ferent methods. In the first one we apply the formula of 
[Horne & Baliunas (1986)[ which resulted in a PAP well 
below 10”^ for the most significant peak. Por a second 
method, we used a bootstrap resampling technique where 
we randomly redistributed the eigenvalue coefficients but 
retained the same times for the observations (e.g., Kiirster 
et al. 119971 1. A 100,000 resampling trials were performed 
but none had reached a power as high as our most signifi¬ 
cant peak and we estimate the PAP to be lower than 10“^. 


The sine-curve periodicity in Pig. 0 is comparable 
to the trend of the Call H&K S-index in Eigure 1 of 
[Metcalfe et al. (2013)[ We see in our data and in the S- 
index measurements that the activity of eEri is higher in 
2009/2010 and lower in 2008/2009 and 2011/2012. The 


rise of activity in the year 2012/2013 is suggestively visi¬ 
ble by [Metcalfe et al. (2013)[ Also the higher variances in 
the years 2011/2012 and 2012/2013 are visible in the Call 
H&K measurements as well. Good agreement is also seen 


for the high activity at the end of the season 2009/2010. The 
overall trend of our results also fits with the S-index mea¬ 
surements from [Jeffers et al. (2014)[ shown in their Eigure 2. 
The behavior of the data from 2009/2010 to 2012/2013 is 
similar to our results of the {B) measurements. The last 
season (2013/2014) seems to be less active than the sea¬ 
sons before. The trend of decreasing activity in 2013/2014 
is also in agreement with the S-index measurements of 
[Jeffers et al. (2014)[ 

5 Conclusions 

Using PCA on a time series of 338 spectra, we were able to 
quantify the magnetic activity variations on e Eri. We com¬ 
puted an average surface-mean magnetic field of 186±47 G. 
The analysis showed variations of the mean surface mag¬ 
netic field on shorter time scales as the rotation period, 
which could be interpreted as spatial inhomogeneities in the 
magnetic field. The lack of a long-term persistent pattern on 
time scales of the rotation period suggests additional tem¬ 
poral inhomogeneities in the magnetic activity of e Eri. 

We confirm the qualitative agreement with the Call 
H&K S-index analysis done by [Metcalfe et al. (2013)[ and 
[Jeffers et al. (2014)| The trend of the Ca II H&K measure¬ 
ments are comparable to our {B) measurements. Our anal¬ 
ysis shows a time variation of the magnetic field that 
is comparable to the shorter magnetic activity cycle es¬ 
timate by [Metcalfe et al. (2013)[ We suggest a period of 
3.0±1.3yr comparable to the 2.95±0.03years determined 
by [Metcalfe et al. (2013)[ By using longer time series spec¬ 
tra of e Eri the method would principally be able to capture 
the longer period as well, if existent. 

The analysis of time series by using PCA is a very 
successful tool to detect variations on cool, slowly rotat¬ 
ing stars. The velocity variations can be separated from 
the magnetic field variations for a detailed analysis of time 
variations of the surface mean magnetic field. The method 
promises to be successful on further stars and will be used 
for analyzing variations of the average surface magnetic 
field on suitable targets in combination with the recently 
presented orthogonal matching-pursuit technique for single 
spectra by Carroll & Strassmeier (12014t . 
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